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Abstract 

^ ' In this paper we study first exit times from a bounded domain of a gradient dynamical system 

\jr^ , Yt = —VU(Yt) perturbed by a small multiplicative Levy noise with heavy tails. A special attention is 

paid to the way the multiplicative noise is introduced. In particular we determine the asymptotics of the 
first exit time of solutions of Ito, Stratonovich and Marcus canonical SDEs. 
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1 Introduction 

In many models of natural phenomena the state of a system is described by a deterministic ordinary differ- 



cntial equation of the form 



r 

Jo 



Y t =y+ \ B{Y s )ds : t>0. (1.1) 



It is often supposed that the vector field B is determined by a function U, so that B = —VU. The function 
U could be called a climatic pseudo-potential [HHHUS] in geosciences, energy potential in physics [TH1HI] or 
profit or cost function in economics and optimization [33] . The potential U is often supposed to have several 
local minima corresponding to the steady states of the system Y. The state space can be decomposed into 
a number of domains of attraction, so that a solution Y t (y) cannot pass from one domain to the other. 

In order to make the models more realistic and allow transitions between the stable states, the system 
(jl.ip is being perturbed by a small random noise, so that (jl.ip turns to a random equation with a small 
parameter. Clearly, the properties of the new random system depend on the interplay between the type of 
the noisy perturbation and the underlying deterministic vector field B. 

Noise can be included into the system in different ways. If the perturbation does not depend on the state 
of the system, one usually speaks about additive noise. If the amplitude of the noise depends on the state 
of the system, one speaks about multiplicative perturbations. 



*Stochastics and Dynamics, Vol. 11, Nos. 2 & 3 (2011) 495-519 



1 



Let for example Z be a regular random process, say, with smooth paths, and F be a smooth bounded 
function. Then the perturbed system with multiplicative noise is described by the random ordinary integral 
equation 



where the last integral is understood in Lebesgue-Stieltjes sense and a positive small parameter e determines 
the noise amplitude. The situation becomes more complicated if one considers irregular perturbations, for 
instance when Z is a Brownian motion. In this case, the differential dZ is usually understood in the sense 
of the stochastic Ito calculus. 

There is a lot of literature devoted to the small noise equation (|1.2[) . both from the point of view of 
Mathematics and applications. The main reference on the large deviations theory and asymptotics of the 
exit times of equation (|1.2p driven by the Brownian motion Z is Freidlin and Wentzell Q2] . In this case, the 
first exit time of X from a domain around the steady state of the underlying deterministic system appears to 
be exponentially large of the order e c / e with the rate C > being interpreted as the energy the Brownian 
particle should have in order to reach the boundary of the domain of attraction. A good exposition of small 
noise properties of Gaussian SDEs with applications can be found in Olivieri and Vares [36] and Schuss |39) . 
Very exact asymptotics of the mean first exit time in the Gaussian case was obtained in Bovier et al. [5] [5] . 

Recently dynamical systems perturbed by small jump noise with heavy tails attracted the attention of 
the physical and mathematical community. The physicists' research focuses on the models incorporating 
a-stable non-Gaussian Levy processes, often referred to as Levy flights. Thus Ditlevsen (UJ [TU] proposed an 
interesting conjecture about the a-stable noise signal in the Greenland ice-core data (see also Hein et al. 
[17| on the statistical treatment of this time series). An enhanced, certainly non-exhaustive list of physical 
references on the first exit problem of Levy-driven SDEs with stable noises includes Chechkin et al. [8] 
and Dybiec et al. QUE]. 

The mathematical theory of large deviations for general Markov processes can be found in Wentzell [47] . 
To our knowledge, in |47] and in Godovanchuk j!6| the asymptotic behaviour of the dynamical systems with 
heavy power tails was considered for the first time. Opposite to the Gaussian case, the behaviour of such 
systems is mainly governed by big jumps. Thus, the exit from the domain occurs with the help of an only 
big jump, and the mean exit time does not depend on the energy landscape of the underlying dynamical 
system, but rather on the geometric layout of the stable states and domains of attraction. Fine small noise 
asymptotics of the SDEs with additive heavy tail Levy noise and their metastablc behaviour was studied 
in Imkcllcr and Pavlyukcvich |23| . The case of light, sub- and super-exponential jumps was considered in 
Imkeller et al. [35]. In his very recent work, Hogele Q2] studied the first exit problem and metastability 
properties of solutions of the infinite-dimensional stochastic Chafee-Infante equation driven by small heavy 
tail Levy noise. 

Coming back to the equation (QHJ) with multiplicative noise, it is necessary to note that the stochas- 
tic integral w.r.t. Z allows interpretations different from the Ito definition, in particular one can consider 
Stratonovich integrals often denoted by o dZ. Even for continuous integrators Z, an interesting question 
arises, namely, which integral fits a specific real world phenomenon, see Arnold [3], Turelli [44], van Kampcn 
[46] . Sethi and Lehoczky [40], Smythe et al. [41] and Sokolov [42] for discussion. Roughly speaking, Ito SDEs 
appear naturally as a continuous approximation of a discrete system, for instance in financial mathematics 
or biology. Due to their nice mathematical properties, they are also the most popular tools in analysis. On 
the other hand, Stratonovich SDEs w.r.t. continuous integrators Z arise naturally as a mathematical ideal- 
ization of dynamical systems perturbed by regular stochastic processes, which takes place in engineering and 
physical sciences. Moreover, Stratonovich integrals enjoy a conventional Newton-Leibniz change of variables 
formula; they are also indispensable for constructing SDEs on manifolds. 

If the integrator Z is a jump process, for instance an a-stable Levy process, the simple Newton-Leibniz 
change of variables formula does not hold any longer even for the Stratonovich integral. To correct this 
situation, the so called canonical SDEs were introduced be S. I. Marcus in [521 133j . 

In this paper we study multidimensional SDEs of the type (jl.2l) . The random process Z is supposed to 
be a multivariate Levy noise with regularly varying (heavy) tails, and the stochastic differential equation 
will be understood in the senses of Ito, Stratonovich and Marcus. Our study treats the first exit time of 
the perturbed system from a bounded domain around a stable attractor of the underlying deterministic 
dynamical system Y in the limit of small noise. 




(1.2) 
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2 Object of study 



2.1 The underlying dynamical system 

Wc start with a n-dimcnsional gradient system generated by a vector field — V£/, 

Y t =y- f VC/(F s )ds, t>0. 
Jo 

We assume that the potential £/ is a C 2 (R n ,R) function with globally Lipschitz continuous first derivatives 

diU(x), 1 < i < ra, and bounded second derivatives didjU(x), 1 s$ i,j < n. We also assume that the 

potential U has a unique global minimum at the origin, U(0) = 0, that is VU(0) = and the Hesse matrix 

(8i8jU(0))2j =1 is positive definite. 

Let Q a R™ be a bounded domain with piece- wise smooth boundary dQ such that e Q. Assume that 

(n(y), — Vt/(y)) < — <5 for y e dQ and some (5 > 0, where n(y) is a unit outward normal at y e cy- 
linder these assumptions is the unique asymptotically stable attractor of the dynamical system Y t (y), 

Yt(y) — » 0, t — > oo; for all y 6 £/ the trajectories Yt(y) do not leave the domain Q. 

Finally let F(x) = (Fij{x))™jZi, x £ R™, be a n x m matrix of smooth bounded real functions with 

Lipschitz continuous bounded derivatives. Let e > be a small parameter. 

2.2 The driving Levy process 

On a filtered probability space (O, J 7 , P) satisfying the usual hypothesis we consider an m-dimensional Levy 
process Z = (Z , . . . , Z m ) with the characteristic triplet (A, v, /i) with a non- negative definite rax m matrix 
A, a vector fi e M m , and a Levy measure with i/({0}) = and ^(1 a \\y\\ 2 ) v{dy) < oo. In other words the 
characteristic function of Z is given by the Levy-Khintchinc formula 

Ee i(\,z t ) = exp ^ _ t ( Ax ^ x ) + f J( e KA,y> _ 1 _ y>l (|s|<1} ) !/(dy)) , A e M m . 
There is also a canonical Levy-Ito representation of Z as a sum 

Z t = nt + W t + z(N(ds, dz) - dsv(dy)) + zN(ds, dz), 

J(0,t] Jo<||z||<l J(0,t] J|ssJ>l 

with W being a Brownian motion with the covariance matrix A and N being a Poisson random measure 
with the intensity measure v. 

To specify the heavy tail property of Z we assume that v is a regularly varying jump measure at oo. Let 
H(u) denote its tail, 

H{u) := u({z e R m : \\z\\ > u}). 

Then for any a > the measure v enjoys the following scaling property: there is a non-zero Radon measure 
m on B(M. \{0}) with m(K \M m ) = so that for any Borel set A bounded away from the origin, ^ A, 
with m(dA) = the relation 

, . viauA) 1 ,. u(uA) 1 . .. 

m(aA) = lim ; / = — hm ; / = — m(A) 
y ' u^+co H(u) a r u^+oo H(u) a r 

holds for some r > 0. In particular, H(u) is regularly varying at infinity with index — r, that is H(u) = 
u~ r l(u) for some positive slowly varying function I. The homogeneity property of the limit measure m 
implies that m assigns no mass to spheres centred at the origin on R m and has no atoms. 

For more information on multivariate heavy tails and regular variation we refer the reader to Resnick 
p] and Hult and Lindskog [T91I20]. 

2.3 SDE with multiplicative noise 

In this section wc briefly remind the main properties of the Ito, Stratonovich and Marcus (canonical) SDEs. 
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2.3.1 Ito SDE 

For simplicity we start in the one-dimensional setting. Let g t be a cadlag adapted stochastic process. Then 
its left-continuous modification g t - is predictable and can be approximated w.r.t. a u.c.p. topology by simple 
predictable processes g^ of the form 



fJi 



(k) 



ffoI{o}(*) + Yi 9j l (T 3 ,T J+1 \{t), 



where = To < ■ • • < Tfc are stopping times and gj are F Tj measurable and bounded. For a Levy process 
Z (or even a semimartingale), the Ito stochastic integral of g w.r.t. Z is then defined as a limit 



J.t k 
9s- dZ s : = lim V g 3 (Z (k) - Z m 
~T 3 i- 1 



in the sense of the u.c.p. topology, see Chapter II in Protter |37j . 

In particular, one can approximate the Ito integral by non- anticipating Riemannian sums. Indeed, con- 



x 



sider a sequence of random partitions = {0 = Tq t[ ^ • • • < < °°} with limsup„ rj: 
a.s., and \\t^\\ := sup fc |rf n) - t^\\ -> a.s. Then 

r-t k 

g s - dZ s = lim V g T (Z w - Z (fc) ) 
Jn fc-*oo r - i j AI r i-i AI 

in the sense of the u.c.p. topology, see Theorem 11.21 in Protter [37]. We refer the reader to Applcbaum [T] 
and Kunita for the theory of stochastic integration w.r.t. Levy processes, and also to Protter [37J for 
the general semimartingale theory. 

Now we introduce the Ito stochastic differential equation with small multiplicative noise. The matrix 
valued function F given, we perturb the equation to obtain 

X t = x- f VU{X S ) ds + e\ F(X S -) dZ s . (2.1) 
Jo Jo 

In the coordinate form this equation reads 

J-t m ft 

d l U{X s ) ds + sY Fij (X S J) dZi, UK n. 

In particular, under above conditions, there exists a strong solution to the equation (|2.1[) , which is a cadlag 
semimartingale and a strong Markov process, [Tl l29l [37] . 

Let / e C 2 (M™, R) and X be a the solution of (|2.1[) . Then the following change of variables formula (Ito's 
formula) holds (Theorem 11.33 in Protter [37]): 



n r-t i " f' 

/(X t ) = f(x) + J] 3 < /(JT i _)dXi + - J 



+ 2 (/(X s ) - /(X s _) - ^ dif{X a -)AX\ 

S<* J=l 

with [X l ,X- 7 ] c being the path-by-path continuous part of the quadratic covariation of X 1 and X 5 . 
2.3.2 Stratonovich SDE 

Let again gt be a cadlag adapted stochastic process and Z be a Levy process, such that the quadratic 
covariation [g, Z] exists. The Stratonovich integral of gt- w.r.t. Z is defined with the help of the ltd integral 

as 



f 9s -odZ s = f g.-dZ a + \\g,Z\\ 
Jo Jo 1 
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The Stratonovich integral can be also interpreted as a limit of Riemannian sums. Let g and Z have no jumps 
in common, that is X!s<t Ag s AZ s = for all t ^ 0, then for any sequence of random partitions t^ 1 ' we have 

n g ( n ) + g („) 
lim ^ ^ 9 Tj " 1 (^ T w Af - Z T c»> At ) = 5s- ° ^ s 

in the u.c.p. topology (Theorem V.26 in Protter [57]). 

The Stratonovich SDE we arc interested in is then written in matrix form as 



X° = x- 



f V£/(X S °) + e f F(X°_) o (2.2) 
Jo Jo 



or in the coordinate form 



(X°)l =Xi-\ diU{X°) ds + sY FiA X °s-) °dZ>, 1 *S i < n. 

Jo 1 = 1 Jo 



It corresponds to the ltd SDE 



X° t =x-\ VU(X° 8 ) ds + e\ F(X° s _)dZ s + *- f F'(X°_)F(X°_) d[Z, Zf s (2.3) 
Jo Jo z Jo 



which in turn should be understood as 

r-t m ~t 

(X°)\ =Xi-\ diU(X°) ds + e^ Fij(X° s _)dZ: 

JO 3 = 1 J° 

J2, ri m 



11 III pt 

£ ~r 2 J dtF^X^FuiX^) d[Z\ Z k f s , 1 < i < n. 

f—1 .• ;„_i JO 



2=ii,fc=i 

Again, from the theory of Ito SDEs we can conclude that the equation (|2.2j) also has a strong solution which 
is a cadlag scmimartingale and a strong Markov process (Theorem V.22 in Protter [57]). 

Stratonovich integrals w.r.t. martingales are generally not martingales. However, when the integrator 
is a continuous scmimartingale, in our case when Z = W + [it, the Stratonovich integral enjoys especially 
nice properties. In this case the solution of (|2.2[) can be obtained with help of the so-called Wong-Zakai 
approximations. Consider the polygonal approximation Z^ n > of the continuous process Z 

Z&> = Z k + (Zu+i-Zh), - < t < k > 0, (2.4) 

n \ " "/ n n 

and a path-wise ordinary differential equation 

A t (,l) =x-\ VU(XW) ds + e f F(X^)dZW (2.5) 
Jo Jo 

with the last integral understood in the Lebesgue-Stieltjes sense. Then in the limit n — > oo, X^ converges 
to the solution X° of (|2.2[) . see Twardowska [55] for a review on the subject and the collection of results. 

Since the processes Z^ are continuous, the approximation (|2.5[) and thus the limiting Stratonovich SDE 
(|2.2[) are often chosen for the description of the real world physical processes. 

Another remarkable feature of the Stratonovich integral consists in a more simple change of variables 
formula. Let / e C 2 (R",R) and X° be an n-dimensional semimartingale. Then (Theorem V.21, [57] ) 



n r-t 

f(X°) = f(X°) + 2 ^/(*°_) o d{Xy 

7 = 1 JO 



+ 2 (/PC) - /(*._) - 2 3 3 -/(jr a °_)A(jr s o y 

»<i j = l 

so that if the pure jump part of 1° vanishes, (A°) d = 0, one obtains the Newton-Leibniz chain rule of 
Stratonovich integrals. In this case, one can construct SDEs on manifolds. 
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2.3.3 Canonical Marcus SDE 



Canonical SDEs were introduced by S. I. Marcus in [3H [33J in order to preserve the flow property and a 
conventional Newton-Leibniz rule for the solutions of SDEs driven by semimartingales with jumps. We start 
with the formulation of the Marcus canonical equation. 

The matrix valued function F given, for any z e IR m we consider the ordinary differential equation 

Ay(V) = F(y(u))z, u =s 0, , 2 g . 

y(0) = x e K n . 

Since all F« are Lipschitz continuous, the vector field F(-)z is complete, that is the solution of (|2.6[) exists 
and is unique for all x e MP and u > 0. Since Fj E C 1 (R",M), the vector field F(-)z generates the flow of 
diffeomorphisms 

<PZ( X ) = y(u,x;z), u^O. 

Wc denote ip z (x) := <pf(x). 

The canonical Marcus SDE is then formally written as 

X°(t) = x - f Vf7(X|) ds + e\ F(X*_) o dZ s (2.7) 
Jo Jo 

where o dZ denotes the Marcus canonical integral. This equation is understood in the following sense: 

X°(t) = x - f VC/(X S °) ds + e\ F(X*_) o dZ s 
Jo Jo 

= x-\ VU(X:) ds + e\ F{X« S _) o dZ% + e f F(X^) dZ A a 

Jo JO JO 

+ ^ (V AZs PC-) - ^ s °- - F(x:_)eAZ s ) 

= x - f Vt/pC) + e f F(X»_) dZ s + ^ f F'(X*_)F(X*_) d[Z, Zf s 

JO JO 1 JO 

+ ^ U az ° (x:_) x:_ f(x:_)saz s 



where the formula after the second equality sign represents the canonical equation in terms of the Stratonovich 
integral, whereas the formula after the third equality sign gives the Ito interpretation. 

Under above conditions, the canonical equation has a unique global solution which is a cadlag semi- 
martingale, see Theorem 3.2 in Kurtz el al. [3U]. Moreover, this solution is strong Markov (Theorem 5.1 in 

TO- 

The jumps of X° occur only when the jumps of Z occur. If Z s does not have a jump at s, AZ S = 0, then 
the trajectory X" moves continuously like the solution of the Stratonovich SDE driven by Z c s . If Z s has a 
jump AZ S at time s, then the trajectory of the solution jumps from the point X°_ to ip eAZs (X*_). That is, 
it flies from the point X°_ along the integral curve of the vector field F(X^_)eAZ s with infinite speed and 
lands at Lp eAZs (X°_). Then the similar movement repeats inductively. It is clear, that if Z d = 0, then the 
Marcus SDE coincides with the Stratonovich SDE. If the noise is additive, i.e. F = const, all three SDEs 
coincide. 

It is necessary to note that the Marcus canonical integral w.r.t. Z appearing in the equation (|2.7[) is not a 
proper integral since it can be defined only for a special class of integrands depending on the driving process 
Z or on the solutions X° of the SDE, namely for processes of the type g(Z t ) or g(X^_), g being a smooth 
function. We refer the reader to Chapter 4.4.5 in Applcbaum [1] and Definition 4.1 in Kurtz et al. [30] for 
details. 

The Wong-Zakai scheme (|2.5|) can be considered also for jump processes Z. In this case, the solutions 
driven by polygonal approximations (|2.4|) converge to the solution of the canonical equation X° in the sense 
of weak convergence of finite dimensional distributions, see Corollary on p. 329 in Kunita |28j . 
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The change of variables formula for solutions of Marcus SDEs has the form of the conventional Newton- 
Leibniz rule. For / e C 2 (R",IR) we have (see Proposition 4.2 in [50] ) 

n r.t 

f(Xf) = f(x) + J] d Xi f(Xt_) o (Xty. (2.8) 

3 = 1 J ° 

Similar to continuous Stratonovich SDEs, Marcus canonical SDEs can be considered on smooth manifolds, 

see [amEojEn]. 

3 The main result and examples 
3.1 The first exit time 

Consider the first exit times of the processes X, X° and X° from the domain Q: 

r x (e) = inf{f > : X t (x) iG], 
r°(e)=M{t>0:X?(x)tg}, 
rt(e) = inf{t > : X?(x) $ Q). 

The main result of the paper is presented in the following theorem: 

Theorem 3.1 Define the sets 

E = E° := {z e R m : F(0)z £ Q} and E° := {z e K m : ip z (0)^Q} 

and suppose that m(E) = m(E°) > and m(£ , °) > 0. Then for any u > —1 and x e Q the following limits 
hold: 

\\ m -Ee- um{E)H{£ ~ 1)T * {£) = limEe- um(B ° )ff(e ~ 1)r - (£) 

= l[ m E,e- um{EO)H{e ~ 1)T * {£) = 1 

e->0 1 + U 

Moreover, there is 7 > such that this convergence is uniform over all x e Q with dist(x, 8Q) > e 7 . 

In other words, the appropriately normalised first exit times converge in law to the standard exponential 
distribution; there is also convergence of all moments. 

Example 3.1 (The sets E = E° and E° are different) Consider a one-dimensional dynamical system 
Y perturbed by a bivariate Levy process Z = (Z 1 , Z 2 ). Let F = i*2) with 

= 7 — , A 2 i 1 and F *( x ) = 7 1\2 , 1 ■ 

(x + iy + 1 (x — i)^ + 1 

In this case according to (|3.1[) . the set 22 = E° is a union of two half-plains (Fig. [U (1.)), 

E = {zeR*: 3* 1 or < -1}. 

1 2 2 1 

The set -E 10 also consists of two halves (Fig. [T] (r.)), 

r = kl 2 : f 1 -or f" 1 *i --e(0,l] 

I Jo ziFi(u) + z 2 F 2 (u) Jo ziFi(u) + z 2 F 2 (u) 

For example, for a bivariate isometric Cauchy process Z with the jump measure v(dz) = ||z||~ 3 cte, z ^ 0, we 
obtain i^e" 1 ) = 2tt£, m(£) * 0.49, m( J B <> ) « 0.45 and hence Et x (e) as Er°(e) * 0.33/e, Er^e) * 0.35/e. 
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Figure 1: The sets E = E° and E 



Example 3.2 (Reduction to additive noise for n = m = 1) In the case n = m = 1, the exit time r° 
can be obtained with help of the change of variables formula (|2.8|) and a trick which was used by Nourdin 
and Simon in |35j . Consider a one-dimensional canonical Marcus equation 

XI = X - f C/'(A S °) ds + e f F{X°_) o dZ s 
Jo Jo 

driven by a univariate Levy process Z and let Q = (—a, b), a,b > 0. Assume that the perturbation is 
uniformly elliptic in Q, that is F(x) > 0, x e [—a, 6], and introduce the function 

Jo ny) F{x) 

Applying the change of variables formula (|2.8p yields the following SDE for the process Y t ° = f(X£): 

Y? = f(Xf) = f(x) f f'{Xt)U'{Xt) ds + e f f(Xt_)F(Xt_) o dZ s 
Jo Jo 

= .f(x)- f f{Xt)U'{Xt)ds + eZ t 
Jo 

Jo 

with B°(y) := — ^ o and y° = /(x). Note that since F is strictly positive and / is monotone 

increasing with /(0) = 0, we can introduce the new effective potential C/ (y) = — J B°(v) dv which is a one- 
well potential with the global minimum at the origin. Thus the process Y° satisfies the SDE with additive 
small noise which has been studied in Imkeller and Pavlyukevich [21] [23] and Imkeller et al. [25] . It is clear 
that A t ° $ (-a, b) if and only if Y t ° $ (/ _1 ( a ). 

For instance, if Z is a symmetric a-stable Levy process with a e (0, 2) and the jump measure v(dz) = 
z|~ 1-Q dz, z 0, then for the first exit time of X° from (—a, 6) we immediately obtain the asymptotics 

limEe- uM ^ £ ° T - = —J— , u > -1, 

e|0 1 + U 

with 

Ar=f I + 1 r 1 . 



4 First exit time of the Ito SDE with multiplicative noise 
4.1 Big and small jumps of Z 

For p e (0, 1) and e sg 1 let us distinguish the small and big jumps of the driving process Z and decompose 
it into a sum 

Z t = L t + % 

with 

th := 2 AZ s l{\\AZ s \\ >e 



'9 

being a compound Poisson process with the characteristic exponent 



Ee i<A,^> = exp ( t f ( e «\y) _ i) v (dy)\ . 



The Levy process L is a process with bounded jumps, ||AL S || ^ e~ p , and thus possesses all moments. 
Moreover, it is a sum of its continuous component L c t = Wt + pt being the Brownian motion with drift, and 
a pure jump part L d . 

Denote by = tq < t\ < T2 < ■ ■ ■ the successive jump times of rj and by </& the respective jump sizes. 
The inter-jump times Tfc = t^. — Tk-i are iid exponentially distributed random variables with the mean value 



ET fc = s- ■■=[ v[dy) = — — -oo, e - 0, 



H{s-p) 



and the probability distribution function P(Tfc ^ u) = 1 — e u/3 ' , u > 0. The probability law of Jk is also 
known explicitly in terms of the Levy measure v: 

P(J fc e A) = p-^iAn {z : \\z\\ > £" p }), A e S(R m ). (4.1) 
4.2 Perturbations by the process eh 

Lemma 4.1 Let p e (0, 1), p £ := ELi and T e := e -9 /or some 9 > 0. There exist Sq = 6o(p) > and 
#o = #o(/°) > such that for all 8 e (0,<5o), 9 e (0, 9q) there are po = po(S) and Eq = eo(p,8) such that the 
estimates 

\eLrJ = 4Ve\\T £ < e 2S and PfleL]^ > e s ) «S cxp(-£- p ) 
ZioZd /or a/Z p e (0,po) a 71 ^ < £ =g £o- 

Proof: Lct us represent the process L as 

L t = L t + p, £ t 

with L being a zero mean Levy martingale with bounded jumps. 
Step 1. We have the following estimates for the mean value p e : 



/4 := ~EL\ = Zi v{dz), 1 < i «S m, 

f ||z|| 2 t/(<fe) = - f w 2 dff(u) s= £~ 2p ff (1), 

Ji<||z||«;e-p Ji 



||Me|| < V^0)^"- 

Consequently, for any p e (0, 1), 9q : = (1 — p)/3 and Jo : = (1 — p)/4 we obtain 

e||/i E |T £ <£ 25 
for all < S < So and < 9 < Oq and e small enough. 
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Step 2. The of quadratic variation process [eL] d is a Levy subordinator 

[^ = £ 2 2l|AL||2 = £ 2 f f ||z|| 2 7V(dz,d S ). 
ss:t ^° Jo<||z||<e-p 

Since the jumps of [eL] d are bounded, its Laplace transform is well-defined for all A e R and equals 

Ee A[eL]? = cxp L f ( e Ae 2 |M| 2 _ ^ 

^ Jo<||z||sS£-'' 

= exp ( - 1 \ (e Ae2 " 2 - 1) dH(uj) . 

For any A > 0, the exponential Chebyshev inequality implies that 



(4.2) 



exp - Ae" - T e 



f 

J0<MS 



(e 



- 1 



)dff(«)). 



(4.3) 



For A = A e := e" 2 " 5 with < 5 < S = (1 - p)/A we have max Q ^ u ^ £ - P \ e e 2 u 2 < A e e 2(1 "'' ) J, as e J, 0. With 
help of the elementary inequality e x — 1 < 2x for small positive x the second summand appearing in the 
exponent in r.h.s. of (|4.3[) can be now estimated as 



T, 



f ( e A ^ 2 " 2 - l)dH(u) 

JQ<u 



I +J 



u 2 dH(u) 



< 2T E \ E e' 



\ 

Jo<ms;i 



u 2 rfi? («; 



+ 2T e A £ e 2(1 - p) 



Jl<U^£-P 



rfi7( u ; 



< CT e \ e e 2 + 2H{l)T e \ e e 2{1 ~ p) 



with C = | $ 0<u<1 u 2 di? (u)| > 0. Consequently, for all < 5 < So and < 9 < 9q we see that the exponential 
inequality 

V{[e~L] d Te > e s ) s£ exp ( - X E s s + CT e X e e 2 + 2H{l)T e X e e 2(1 ' p) ) < e"^ 2 



holds for e small enough and the Lemma holds with p e (0, 5/2). ■ 

Lemma 4.2 Lei p e (0, 1) and (<7t)t>0 ^ e a bounded adapted cddlag stochastic process mit values in M m 7 
T £ = e -61 , > 0. There are Sq = S (p) > and 9 = (p) > such that for all S e (0, do) and 9 e (0,# ) 
£/iere are po = Po{p,S) and Sq = £q(S) such that the exponential estimate 



P ( sup e 



.,=1 J o y 



/io/ds /or cdZ p e (0,po) ari d < e < £o- 

Proof: Step 1. Suppose that sup 4;s0 ||g t || C for some C > 0. Consider the one-dimensional martingale 



J=l J ° 



dLl. 



By construction |AM t | < Ce~ p . We estimate the probability of a deviation of the size e s of eM t from zero 
with help of the exponential inequality for martingales, see Theorem 26. 17(i) in Kallenberg [26]. Indeed for 
any 8 > and 9 > we have 



P( sup I eM t I 3s e 5 ) < P( sup |eM t | 3s £ 5 [eM] Te < 



_4<5 



P([eM] Te > 



.4,5 \ 
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Inspecting the proofs of Lemma 26.19 and Theorem 26. 17(i) in Kallenberg [55] we get that for any A > 



P sup eM t > e e 



[sM] Te < 



-4i> 



e 



-XE S + X' 2 h{XCs 1 ~ p )e 4S 



with h(x) = — (x + ln(l - x) + )x 2 . For any < 8 < 8i := (1 — p)/2 we set A = X £ = e 25 , so that 
h^eCe 1 ^ 9 ) — > 1/2 as e — » 0. Hence we obtain the estimate 



P( sup eM t > e £ 

1 ts=T e 



[eM] Te < e 



4(5 



which holds for e small enough and p e (0, 5/2). An analogous inequality holds for inf^^ sM t . 
Step 2. There is a constant Ci > with 

[sM] t = f 5s 2 d[eW] s + f 5s 2 _ *S C\(e 2 t + [eL]*), t > 0. 

Jo Jo 

For sufficiently small e we obtain the estimate 

P([eM] Te > e 45 ) < P([ei]£. > e 55 ) + P(C l£ 2 T e > e 55 ). 

For e small enough the second summand vanishes for all < 5 < 82, < 8 < 82 with < 82 + 662 < 2. 
The first summand is bounded by e~ £ from above due to Lemma 14.11 for < 5 < 83 some < 8 < 83 
and < p < pi(p, 8) and e small enough. This the statement of the Lemma holds with (5o = min{<5i, 82, 83}, 
6*o = minj^i, 82, 63} and po = mm{S/2,pi(S)}. ■ 

Lemma 4.3 For any p e (0, 1) there are 70 > and p > such that for all < 7 70 



xeg 



supP sup \\X(t,x) - Y(t,x)\\ Ss e 7 «S e 



o^t<Ti 



Proof: Step 1. By assumptions on the potential U, any deterministic trajectory Y t (y), y e Q, reaches a 
small fixed neighbourhood of the origin in some finite time. After entering this small neighbourhood, the 
trajectory Y t (y) is attracted to the origin with the speed approximately proportional to Ci||Yt||, Ci > being 
the smallest eigenvalue of the matrix -5-^= — U(x)\ n . This allows us to estimate of increase rate of the time 

dxidxj \ '\x=0 

a trajectory Y t (y), y e Q, needs to reach some e -neighbourhood of the origin. Indeed, for any 6 > the 
following inequality holds true for any O^e^eo, y £ G and £0 > small enough: 

28 



\\Y t {y)\\ < e A , t>V e 



C x 



In; 



(4.4) 



Step 2. Here we show that on the time intervals up to V £ , the random trajectory X s (x) does not deviate 
much from the deterministic solution Yt (x) with the same initial value in the absence of big jumps of the 
driving process Z. 

For x G Q, with help of Gronwall's lemma we estimate 



\\XtAVeATi_- — YtAVcATi-l 

rt a V e aTi - 



Jo 



X.-YJ ds + e 



Jo 



F{X 8 _)dZ s 



sup \\X t -Y t \\ <e CLi ^ sup e f F(X s _)dL 6 

0sSt<V e ATi 0^t^V e Jo 



Recalling the definition of V £ in (|4.4[) and taking into account that V e < T e 
any > we get for any 8 > that 



sup 

0^t<V e aTi 



A t -y t || >s° <P sup e 



Jo 



5* £ 



for e small enough and 

) 



(1+^)4 



<P( sup el f F(X s _)dL s + C 2 eT e \\ij, e \\ ^ e {1+ ^r )6 ) 

^ Osgt=ST E I Jo ' 
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with some C2 > 0. With help of Lemmas 14. II and 14.21 we hnd that there are > and 9± > such that the 
last probability is smaller than cxp(— e~ p ) for all 6 e (0, S±), 9 e (0, 9\) and p e (0,pi(<5)). 

Step 3. In this Step we exploit the attractor property of the origin and show that in the absence of big 
jumps of the driving process Z the random path Xt(x) with the initial value close to the origin does not 
deviate much on the polynomially long time intervals T e . 

Consider the function f(x) = ln(l + U(x)) > 0. For ||x|| small, one can estimate ci||cc|| 2 ^ f(x) ^ C2 1| re || 2 
for some positive C\ and c 2 . Furthermore, the derivatives dif(x) = j^jj^j, Sidjf(x) = 2 s * u ( x )^u(x)-d^djU(x) 
are bounded due to assumptions on U. 

We apply the Ito formula to the process f{X t ): 



0^/(1^^-) 

n r t a T E a Ti - 



m + S 

•i=l J 



d l f{X s _)dXl + 



n r t a T e a Ti - 



1 r 

2 U=i Jo 



d^/pC.-) JC J '] 



/(*) - f 
Jc 



+ 2 (f(x s ) f{x a j) 2 aj(x s _)Ax: 

s<(aT b aTi i=l 

tAT » AT i-||VC/(X s _)|| 2 



£ 



l + f/(X s _) 

2 ra m r t a T E a 7\ - 



1 + C/(X S _) 



U III r t, 

■s s f 

i,j = l fc,{=l J0 



d i d j f{X.-.)F ik { L X s J)F jl {X s _) d[Z k , z 1 ^ 



+ 2 (/(^)-/(^-)-£^/(X s _)AJf: 



s<t aT £ aTi 



We note that the first integral in the last formula is non-negative, the integrands in the Ito integral w.r.t. Z 



and in the integrals w.r.t. [Z k , Z l ] c are bounded, the quadratic covariations satisfy [Z k , Z 1 ]^ = [W k , W 1 ^ 
<t kit, and finally the estimate 

n 

2 /pr.) - /(x s _) - 2 dif{x a J)6JC t 



, n f l 



(1 - v)d i d j f(X s _ + vAX s ) dv 



<c 3 2l|AJr 8 | 2 = C3[Jf] t d 



holds with some C3 > 0. Furthermore, since F is bounded the estimate 

[X] d t < C 4 [Z]f = C 4 [L]f 

holds for some constant C4 > and < t < T\. Combining these estimates and denoting g(x) = V f^^fr^ 
we obtain the following estimate with some positive constant C5, £ small enough and ||x|| =S e s , S > 0: 



0</(X tA T E AT 1 -)*£C 5 (£ 25 + £ Sup 

Os£tsST. 



f 3(^ s 
Jo 



+ £||^||r e + £ 2 T £ + £ 



Let the estimates of the Lemmas 14.11 and 14 . 2 1 hold simultaneously for S e (0, £2), 9 £ (0, 6*2) and p e (0, P2(5)) 
for some positive £2, 92, P2($) and £ small. Then for S e (0,(52/3), 9 e (0,#i) and £ small we get 



sup |jJT t || 

0sS4<T s aTi 



s£ P ( £ 2<5 > £ 3<5 



+ P 



P [e\\^\\T £ >e 



sup 

0^t^T e 



P(e*T e > £ 35 ) + P( > e M ) *S e 



f g{X s ^)dL a 
Jo 



> £ 
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forpe (0,p 2 (<5)/2). 

Step 4. Combining the estimates of Steps 1, 2 and 3, we extend the estimate of the Step 3 to all initial 
values x e Q: 

sup \\X t -Y t \\^e s ) 



^OsSt<T e ATi ' 

for S e (0,^3), $ e (0, #3), p e (0,P3(S)) and e small enough. Here, ^3 = min{<5i, 62/%}, 6*3 = min{0i,#2} and 
p 3 = mm{p 1 (S),p 2 (S)/2}. 

Step 5. In this final Step we extend the estimate of the Step 4 from the time interval [0,T e a T\) to the 
time interval [0, T\). 

Denote X L the solution of the SDE (|2.ip driven by the process L. Clearly, X t = X^ on the event 
{t < Ti}. Let x e Q and fc > 1, then for any <5 > and 9 > we have 

P( sup pT t -y t || >£ 5 ) ^Pf sup \\X t -Y t \\ ^s s ) +P(T 1 > kT E ) 

v te[0,Ti) ' v te[0,feT e AT!) 7 

sSPf sup ||X t L -y t || +P(Ti >kT e ). 

Moreover for e small enough we have ||y(T e ,x)|| < e 2d for xe^. Denote 

A 3 = { sup IX^ - Y (t - jT s ; Xf Tt ) \\ < e s }, < j < fc - 1. 

te[jT £ ,(i+l)T £ ] 

In particular, the probability of Aq = {sup te [ T i \\X^ — Y t \\ > e s } was estimated in Step 4. Further, for any 
k>l 

k-l 

S 



f|As{ sup ||X t L -r t || <2e 



J=0 - te[0,fcT eJ 

Consequently 

fe— 1 

sup |X t L -y t |>2 e i ) <P(|J 



' (Al u (A Ai) u • • • u (A • • • Au^A^S) (4.5) 



fc-i 

P(^, ||Xf T J 6 5) < A 

-2r 



< 2 P(A|, |LY/ T J E 5) *S fcsupP(^). 

j=0 



For k = k e = [s ] and any 9 > we have 

P(Ti Ss fc e T e ) = e - fe » T -A ^ exp(-e rp - e - 2r /(e- p )) < e^ P . 
for all < p < P4 := (2 — p)r and e small. On the other hand, (|4.5[) and Step 4 yield 



P( sup \\Xf -Y t \\> 2e 6 ) e-^ r e- £ e 

te[0,fcT E ] 

with any < p < pz{5). Consequently, the statement of the Lemma holds for any 0<7<<53,0<p< 
min{p3 (7)72,^)4} and e small enough. ■ 



4.3 The first exit time of solutions of the Ito SDE 

Having established the key estimates about deviations of the random trajectory Xt(x) from the deterministic 
path Yt(x) on random time intervals between big jumps of the driving process Z we can calculate the 
asymptotics of the Laplace transform of the first exit time. The proof here goes along the lines of the 
one-dimensional case considered in Imkcllcr and Pavlyukevich |22| and the multivariate case of a dynamical 
system driven by a multifractal a-stable noise considered in Imkeller et al. [24] . For the sake of completeness 
we briefly sketch the main idea of the proof. 
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The argument is based on the concept of the one big jump which is often used in the study of heavy tail 
phenomena. Roughly speaking, it can be shown that under certain conditions the small perturbations of the 
dynamical system Y due the process sL can be neglected, and the exit from the domain occurs with high 
probability at one of the jump times r k . Just before the time r k the solution X stays in a small neighbourhood 
of the stable point, so the exit occurs if the jump eJ k is large enough, namely if F(X Tk _)eJ k F(0)eJ k £ G- 
The events {eJ, $ E) = {F{Q)eJ x e G},. ■ ■ , {eJ k -i $ E} = {F(0)eJ k ^ e G) and {eJ k e E} = {F(0)eJ k $ G} 
are independent and build up a geometric sequence of events. Their probabilities can be calculated in the 
limit of e — > with help of the scaling property of the jump measure v. 

The statement of the main theorem follows from the small noise estimates from below and above of the 
Laplace transforms of the normalised first exit times. Here we consider the less complicated estimate form 
below. 

For any u > — 1, with help of the formula of the total probability we have 

Ee -um(i?)ff(£-> x (e) > ^ ■E 1 \e- um{E)H ^ 1)Tk l{a = r k } . (4.6) 
fc=i 

For any 5 > small enough denote G~ S := {x e G '■ dist(c5G, x) > 5} the inner part of G and G +S = G S '■ = 
{x e M™ : dist(G,x) <5} be the outer ^-neighbourhood. For k > 1, the strong Markov property allows to 
write 



E 



e - tl m(_E)_f/(e- 1 )r fc 



I{<7 = Tfc}] 

E ^ e -um(E)H(e-i)T kI{Xt e g^ e [ , Tfc ),X Tfc * Q} 
fe-1 

v[Yl e ~ um{E)H{E ~ 1)T3l { x t+T 3 -i eg,te [0,2}]} 

x e-™W H ^ T *I{X t+Tk _ 1 eQ,te [0, T k )}l{X Tk $ G}] 
inf E\e- um ^ H ^ Tl I{X t (y) e g~ s ,te [0,^]}]) ~* 

eQ- s L J/ 

inf Ele-^W^^Xtiy) eQ,te [O,^)}!^ $ g}\ 



(4.7) 



Lct 7 > be such that the estimates of the Section l4~21 hold. We set 5 := 5(e) = e 7 . The exit from the domain 
with a big jump eJ\ occurs when F(Xt 1 -)sJi $ G- Further, sup 0;St<Tl \\X t — Y t \\ sg ig 7 with probability 
exponentially close to 1 (Lemma 14. 3[) , Y t (x) reaches a ^^-neighbourhood of the origin during the relaxation 
time V e = C(|lne|), and 2\ > V e with high probability. Taking into account that Ti is exponentially 
distributed with the parameter /3 e we calculate the Laplace transform of m(E)H(e^ 1 )Ti explicitly, namely 

^ e -u m{ E)H{e-- ) T 1 = Pe = _J_ m[E) Bfll 

f3 e + um{E)H(e- 1 ) 1 + ua E ' v H(^ p ) 

Recalling the probability law of big jumps (|4. 1 [) we see that for e small enough 



whereas for any 6' > 



P(F(0)eJ k i G) = P(eJ k e E) = ^M- (4.8) 

Pe 

a e {l-5')< V -^-^a e {l + 5'). (4.9) 

Pe 

To obtain the final asymptotics we have to estimate carefully the perturbed the exit probabilities P(F(y)eJi e 
G~ el ) and P(F(y)eJi $ G) uniformly over \y\ < e 1 . This is achieved with help of the continuity of the 
function (y, z) i— > F(y)z both in y and z. Indeed, for any 5' > we can choose R > big enough, such that 
the estimate 
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holds for e small. Further, the function F(y)z is uniformly continuous in z in the ball ||z|| ^ R and is 
continuous in y at the origin. Using the scaling property of the jump measure v and the fact that the 
limiting measure m has no atoms we show that uniformly over (|y|| < e 1 

\P(F(y)sJ k i Q ±£ \\\eJ k \\ < R) - P(F(0)eJ k $ G, \\eJ k \\ *g R)\ ^ 5 



4 H(e-py 
and 

P(F(0)sJ k i G) - P(F(0)sJ k i Q, \\eJ k \\ <R)^ jffeA 
Finally for any 5 > we choose 6' > small enough to get the uniform estimates 

inf P^eJ, e g- £l ) > 1 - m (£?)^f_^(i + <5) = 1 - a £ (l + <5), 
inf P(F(y)eJi £ 5) > m(E)^f^-(l - S) = a £ (l - S). 

for e small enough. 

Following the lines of the proof of (22] and [24] , for any 5 > and £ small we can also obtain the 
multiplicative estimates for the Laplace transforms for any u > — 1: 



inf E 



e -um(E)H(e-i)T 1I{Xtiy) £ g-^f e [0 , Tl ]}] > 1 SeH+H 



1 + HO E 



inf V\e-™W H ^^I{XM e £T'"\f e [0,T 1 )}I{X Tl * £}1 ^ * } . 
Summing up the terms from (|4.7|) over fc > 1 yields the estimate 

Ee -™(£)fl( E ->«( E ) ^ Qe(l-^) y / l-a £ (l + ^) ^-i 

1 + ua e V 1 + ua e 
k — i 

5, i _ ^ 



l + u + 5 1 + u 

for some C(it) > and e small. 

The upper bound for the Laplace transform is technically more involved since it additionally demands 
careful estimates of the probability to exit from the domain due to small jumps during the inter-jump intervals 
of the compound Poisson process rj. These estimates are obtained analogously to the one-dimensional 
and multi-dimensional cases studied in |22[ I24j and finally lead to the uniform convergence of the Laplace 
transform over x e g~ £ as e — ► 0. 



5 First exit time of the Stratonovich SDE 

Recalling the Ito form of the Stratonovich SDE (|2.3[) . we reduce the exit problem of X° to the Ito case. 
Indeed, in the argument of the Section 14.31 we have to take into account the Stratonovich correction term 
^- J Q F'(X°)F(X°)d[Z, Z] c which is a Lebesgue integral whose absolute value increases at most as Ce 2 t for 
some C > 0. It is clear that adding this term to the equation does not influence the estimates of the section 
14.21 Thus the result follows immediately, and we obtain the same asymptotics of the first exit time as in the 
Ito case. 



6 First exit time of the Marcus (canonical) SDE 

The analysis of the canonical Marcus SDE can also be reduced to the Ito case. As in Section l4~3l above let 
us distinguish between big and small jumps of Z. Since the processes rj and L = Z — rj are independent, the 
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Marcus equation can be rewritten in the Ito form as 

ft rt 



X° t =x- \ VU(X*) ds + s\ F(X°_) o dZ c s + e\ F(X°_) dL d s 
Jo Jo Jo 

+ 2 (V AL * (xt-) x:_ + f(x:_)salA + 2 U Ai *(x:_) xt 



Let us estimate the small jump correction term in the Marcus equation. The jumps of the process eL d arc 
bounded in absolute value by e 1 ~ p . The mapping u *-* <fi% is C 2 (]R,R n ), so the Talyor expansion yields 

<p z (x) = y(l,x;z) = y(0,x;z) + —y(0,x;z) + ^^vi ^ x '> z ) 
= x + F(x)z + R(x,z), 6(0,1), 

where 



2 yeK" I ^ fez 



1 k sS ra. 



Since all i^j are bounded with bounded derivatives, we obtain the estimate 

\cp z (x) -x-F(x)z\ ^C\\z\\ 2 
with some absolute constant C > 0. This leads to the inequality 

| 2 (> Ai »pC_) + W_)eALf)| < C^e'llALfH 2 = C £ 2 [L]f, (6.1) 

so that this summand is small due to Lemma 14.11 Thus we are again in the setting of the deterministic 
dynamical system Y perturbed by a small noise process 



ft 

e 



\ F{Xt_) o dZt + s f F(X:_) dL d s + ^ (X:_) X:_ + F(X s «_) £ ALf) (6.2) 

JO JO s <; t V 7 

and a big jump process ip eJk (X* _) — I° t _ . 

The arguments of the Section 14.31 for the proof of the Ito case can be applied to the Marcus canonical 
equation. First, due to the estimate (|6.1[) and Lemmas 14.11 and 14.21 we obtain the exponential estimate of 
the Lemma POl for the solutions of the canonical Marcus SDE X°. 

Then we again exploit the concept of the one big jump, that is we show that the exit from the domain 
Q occurs with high probability at one of the jump times Tfe. Just before the time the solution X° stays 
in a small neighbourhood of the stable point, so the exit occurs if the jump eJk is large enough, namely if 
<pt J x(X Tk -) * ^ Jfc (0) i g. The events {e.h $ £°} = {^ Jl (0) e £},..., {eJ k -i i E°} = {^^(O) e'^} 
and {eJfc e E } = {ip sJk (0) $ G} are independent and build up a geometric sequence of events. As in the Ito 
case, for any 6' > their probabilities can be calculated in the limit of e — * as 

Fine estimates for the perturbed the exit probabilities P(<p Ejl (y) e Q^ 6 ' 1 ) and P(tp eJl (y) $ Q) are also 
obtained analogously to the Ito case. Indeed, one can see the mapping (y, z) i— » F(y)z as a particular case 
of the mapping (y,z) <p z (y) appearing in the Marcus equation. Again, ip z (x) is continuous in y at y = 
and is uniformly continuous w.r.t. z in the ball ||z|| < R with some R big enough. Thus the argument of the 
Section l4~3l can be repeated directly with ip z (x) instead of F{y)z and E" instead of E. Consequently for any 
5 > we obtain the uniform estimates 

inf P(^(y) 6 ep") > 1 - m(EO)§f^r(l + S), 

inf P(vJ eJl (y) i G) > m(^)f^(l - 5) 

\\y\\ae—i H(e p) 

for e small and hence also the estimates for the Laplace transform. 
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